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10.3 IF ESTIMATION FOR MULTICOMPONENT SIGNALS 0 

10.3.1 Time-Frequency Peak IF Estimation 

There is a wide range of applications where we encounter signals comprised of M 
components with different IF laws / m (t) and different envelopes in additive 

noise. It is often desired from such an observed signal, to determine the number of 
components M, the IF law of each component and the corresponding envelope a m (t). 
This can be achieved by representing the observed signal z(t) in a time- frequency 
(t-f) domain and use time-frequency filtering methods to recover the individual 
components [1]. Another approach involves extending algorithms for IF estimation 
of monocomponent FM signals to the case of multicomponent signals and design 
an algorithm that simultaneously tracks the various IF components of the observed 
signal [2,3]. Both approaches require the use of time- frequency distributions (TFDs) 
with very specific properties such as high time-frequency localization of the IF 
components and high reduction of cross-terms interferences. 

The basic concept of instantaneous frequency is described in Pt. 1 of reference [4] 
and in Chapter 1 of this book. Methods of IF estimation are reported in Pt. 2 of [4] 
and in Chapter 10. Essential results are reproduced below, for greater clarity. 

10.3.1.1 Spectrogram Peak IF estimation 

Various approaches for IF estimation of monocomponent signals exist [4]. Most 
of these algorithms are suited to a particular class of signals, and both fixed and 
adaptive algorithms have been proposed. Our aim here is to approach the problem 
from a general viewpoint in order to define a general IF methodology that would be 
suitable for the largest class of signals found in practical applications. To illustrate 
this approach, we thus consider from the outset multicomponent signals in additive 
noise, which can be expressed as follows: 

M M 

z(t) = ^2 z ™{t) + e(t) = a m{t)e^ m{t) + e(t) (10.3.1) 

771 = 1 771=1 

where the amplitudes a m (t) are the component amplitudes, are the compo- 

nent phases, and e(t) is a complex-valued white Gaussian noise process of indepen- 
dent and identically distributed (i.i.d.) real and imaginary parts with total variance 
of. The individual IF laws for each component are given by [1]: 

fm{t) = ^ d ^r L = (10.3.2) 

Z7T at 
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Fig. 10.3.1: Left: The spectrogram of a bat signal using a small analysis window. Right: The modified 
B-distribution of the same signal with parameter (3 = 0 . 05 . Total signal length is N = 400 and sampling 
interval At = 1. The spectrogram cannot show the weakest component. 



A conventional approach to represent and analyze such signals for IF estimation 
is to take the spectrogram of z(t) and search for the peaks in the t-f domain (see 
Article 10.1). Curves formed by a continuum of these peaks describe the IF laws of 
the individual components of the observed signal z(t ), as illustrated in Fig. 10.3.1 
using a bat signal. Analytically, this can be expressed as follows: 

fmit) = arg[max p m (t, /)]; 0 < / < f s /2 (10.3.3) 

where p m (t , /) is the m th peak of the spectrogram. 

This spectrogram based approach has several advantages: it is easy to under- 
stand, easy to use, and there are no cross-terms producing unwanted interferences. 
A major disadvantage though is that the time-frequency resolution of the spectro- 
gram for closely spaced components is often poor, especially if one of the compo- 
nents is much weaker, as illustrated in Fig. 10.3.2 for a two-component linear FM 
signal with one component weaker than the other. Fig. 10.3.2 also compares the 
performance of the spectrogram with the modified B-distribution (MBD) that was 
designed specifically for multicomponent IF estimation, as discussed later. 

10.3.1.2 Peaks of WVD, PWVD, and RIDs 

To improve upon the resolution of the spectrogram, various TFDs were proposed 
for IF estimation, one of the most important being the Wigner-Ville distribution 
(WVD). IF estimation using the peak of the Wigner-Ville distribution (WVD) is 
optimal for linear FM signals with high to moderate signal-to-noise ratios (SNRs) [4] , 
but its performance degrades significantly at low SNRs, and in this case the cross 
WVD (XWVD) peak can be used as an IF estimator [5]. For polynomial FM 
signals it was shown that the polynomial WVD (PWVD) gives the best performance, 
especially at high SNRs (see [6] and Article 10.4). However, both WVD and PWVD 
suffer from cross-terms when used to analyze multicomponent signals. These cross- 
terms generate artifacts that obscure the (t-f) representation of the signal, leading to 
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Fig. 10.3.2: Performance comparison between the spectrogram Spect(t, /) and the modified B- 
distribution (MBD) for (3 = 0.06 using a two-component noise-free linear FM signal at the sampling 
instant n = t/ At = 64. Total signal length is N = 128 and the sampling interval is At = 1. The right 
component is five times larger in amplitude than the left component. Left: Spectrogram with small 
analysis window length (A = 23). Right: Spectrogram with large analysis window length (A = 83). 
In both cases the spectrogram fails to resolve the two components. In addition, time resolution is bad 
for a large window length. 



the development of reduced interference distributions (RIDs) to remedy the problem 
[7]. Straightforward IF estimation using the peak of RIDs give an IF estimate that 
is biased from the true IF law, and this bias is different for different RIDs. Although 
reduced, cross-terms still exist and can obscure weak components, hence the need to 
define special purpose RIDs with efficient cross-terms reduction, high time- frequency 
resolution and minimum bias from the true IF laws, such as the MBD [3,8]. 

10.3.2 Properties of IF Estimates Based on Quadratic TFDs 
10.3.2.1 IF Estimates and Window Length 

We consider an analytic signal z(t) of the form z(t) = ae + e (f) where the 
amplitude a is constant, and e(t) is a complex-valued white Gaussian noise with 
independent identically distributed (i.i.d.) real and imaginary parts with total vari- 
ance of . The IF of z(t) is given by eq. (10.3.2), and it is assumed to be an arbitrary, 
smooth and differentiable function of time with bounded derivatives of all orders. 
The general equation for quadratic time- frequency representation of the signal z(t) 
is given by [1] 

Pz{t,f)= T\ G{t,T)*K z (t,T)\ 

( t ) 

where G(t,r ) is the time-lag kernel, K z (t,r) = z(t + | )z*(t — |) is the signal 

kernel or the instantaneous autocorrelation function (IAF), and * denotes time 

(t) 

convolution. For smoothing and localization on the IAF we apply a window function 
Wh(r) = on the instantaneous autocorrelation K z (t,r), where w(t) is a 

real- valued symmetric window with unity length, i.e. , w(t) = 0 for | t |> hence 
the window length is h. 
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The TFD is now dependent on the window length h as follows: 

Pz,h(t,f)= T \w h (T)G(t,T) * K z (t,r)\. (10.3.4) 

T ^f (*) 

If Pz,h(t, /) is discretized over time, lag, and frequency then we have 



Pz,h(n,k) 



N a -1 N a -1 

Wh(rnAt)K z (lAt, 2mAt)G{nAt — lAt , 2 mAt)e~-’ 2n 2 N a 



l=—N s m=—N s 

(10.3.5) 



where 2N S is the number of samples and At is the sampling interval. 

The IF estimate is a solution of the following optimization 

f h (t) = arg[max p z>h (t, /)] ; 0 < / < f s /2 (10.3.6) 



where f s = 1/At is the sampling frequency. 



10.3.2.2 Bias and Variance of the IF Estimate 

By extending the results in [9] , the estimation bias and variance are found to be [3] 
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(10.3.7) 
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(10.3.8) 

where A cf>(t, r) = 4>{t + r/2) — cj>(t — r/2) — r</> (t). 

Equations (10.3.7) and (10.3.8) indicate that the bias and the variance of the 
estimate depend on the lag window length h for any kernel G(t , r). To see how the 
bias and the variance vary with h, asymptotic analysis as At — > 0 is necessary for 
the chosen TFD. 



10.3.2.3 TFD Properties Needed for Multicomponent IF Estimation 

The results above indicate that a general method for IF estimation of multi- 
component FM signals in additive Gaussian noise that is based on quadratic 
time-frequency distributions requires the quadratic TFDs to satisfy the following 
conditions [3]: 
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(1) p(t, /) should have a high time-frequency resolution while suppressing cross- 
terms efficiently so as to give a robust IF estimate for mono- and multicomponent 
FM signals. 

(2) p{t 1 /) should enable amplitude estimation for the individual components of the 
signal, as the amplitude is necessary for evaluating the variance of the IF estimate 
for each component [2,3,9], and to allow for the reconstruction of the individual 
components of the signal. 

(3) the choice of the lag window length should lead to a bias-variance tradeoff (see 
eqs. (10.3.13) and (10.3.14)). 

Although some TFDs, like Choi-Williams distribution (CWD) and the spectro- 
gram, can satisfy some of these conditions, they do not meet the second requirement, 
i.e. allowing direct amplitude estimation. The design of TFDs which satisfies all of 
these required properties is considered next. 

10.3.3 Design of Quadratic TFDs for Multicomponent IF Estimation 

10.3.3.1 Desirable Time-Lag Kernel 

A TFD referred to as the B-distribution (BD) was proposed and shown to be su- 
perior to other fixed-kernel TFDs in terms of cross-terms reduction and resolution 
enhancement [10]. As it does not allow direct component amplitudes estimation [2], 
as per the second condition on TFDs required for multicomponent IF estimation 
listed above, the BD kernel was modified as [3] 

G{t, t) = G f}{t) = fc ( g/cosh 2/3 (t) (10.3.9) 

where (3 is a real positive number and kp = r(2/3)/(2 2/3 - F 2 (/3)), T stands for the 
gamma function. This modified B-distribution MBD(t, /) is also referred to as the 
hyperbolic T-distribution (HTD) in [8] . 

10.3.3.2 Relevant Properties of the Modified B-Distribution (MBD) 

Most of the desirable properties of time-frequency distributions relevant to IF 
estimation (as explained in [1] and [4]) are satisfied by the MBD kernel. In 
particular, realness, time-shift and frequency shift invariance, frequency marginal 
and group delay, and the frequency support properties are satisfied. The time 
support property is not strictly satisfied, but it is approximately true [3]. The three 
required conditions listed in Subsection 10.3.2.3 are discussed in detail below. 

(1) Reduced interference and resolution: This property is satisfied by MBD. 
For example, consider the sum of two complex sinusoidal signals z(t) = Z\(t) + 
Z 2 (t) = aie^ 2n ^ lt+6l ' > + a^e^ 2 ’ 1 '^ 2 * -1 " 6 ' 2 ) where a\, 02 , 0\ and 0 2 are constants. The 
TFD of the signal z(t) is obtained as [3] 

MBD(f, /) = al 6(f - h) + a 2 2 S(f - f 2 ) + 2a 1 a 27/3 (t) S[f - (/1 + f 2 )/ 2] (10.3.10) 
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where 7 p(t) = |T(/3 + jn(fi - / 2 ))| 2 cos (2tt(/i - f 2 )t + 6 >i - 6> 2 )/r 2 (/3). The cross- 
terms are oscillatory in time and depend on the frequency separation between signal 
components. If f\ and / 2 are well separated then the term |T(/3 + jn(fi — / 2 ))| 2 
can be substantially reduced, while T 2 (/3) can be made high if (3 is small. When /1 
and / 2 are not well separated, the MBD still performs better than most quadratic 
TFDs (see Article 7.4). 

(2) Direct amplitude and IF estimation: The MBD allows direct IF estimation 
by peak localization, i.e. , at any time instant t, it has an absolute maximum at 
/ = 27 for linear FM signals . For non-linear FM signals this estimate is biased, 
but this bias can be accounted for in the adaptive IF estimation, as presented next. 

For an FM signal of the form z(t) = a the MBD is approximated by [3] 

MBD(f, /) «| a | 2 J G p{t-u) 6 [^'(u) - f]du = \a\ 2 G 0 (t - ip(f))ip' (f) 

(10.3.11) 

where ip is the inverse of ^(p , i.e., p^cp ( ip(f )) = f. Assuming that xp (/) is not a 
highly peaked function of f and knowing that G 0 (t — ip(f)) is peaked at t = 
the absolute maximum of MBD(t, /) for any time t would be at ip{f) = t, or 
/ = 7p(p (t), which is the IF of the FM signal z(t). For non-linear FM signals, the 
energy peak of the MBD is actually biased from the IF because of the extra term 
Y^P= 3 (k odd) (p ik) {u). The major contribution in this term is due to <p^ 3 \u). 

Therefore at the instants of rapid change in the IF law the bias is not negligible and 
eq. (10.3.11) would not be an accurate approximation to the MBD unless suitable 
windowing in the lag direction is used. 

For linear FM signals we have <p^ k \t) = 0 for k > 3. Assuming that z(t) = 
ae i2ir(/o*+^-« ) ; where f 0 and (3 0 are constants, we have 

MBD(f, /) - jW\ 2 G 0 (t - jU - f 0 j) (10.3.12) 

which has an absolute maximum at / = f 0 +f3 0 t , the IF of the linear FM signal z(t). 
As (3 0 — i - 0, i.e., z(t) approaches a sinusoid, we have MBD(t, /) — ► \a\ 2 S(f — f 0 ), 
in accordance with eq. (10.3.10). 

As for amplitude estimation, eqs. (10.3.11) and (10.3.12) indicate that the MBD 
can support amplitude estimation [3]. 

(3) Asymptotic Formulas Using MBD: The asymptotic formulas for the vari- 
ance and the bias as At — > 0 using a rectangular lag window are given by [3]: 
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and 



where A (f) = f^ 2 \t + n) + /( 2 )(f — n), is the second derivative of the IF, 

and sup|/( 2 )(t)| < Mi- For small h, the optimal window length that minimizes the 

t 

mean squared error is obtained by extending the result in [6] as: 

h r 1800^(1 + ^) 

optW ~ U 2 |a| 2 (/( 2 )(t)*l/cosh 2/3 (t)) 2 
(*) 

Thus, the optimal window length depends on the second derivative of the in- 
stantaneous frequency which is time and signal dependent. Eqs. (10.3.13) 

and (10.3.14) indicate that the variance and bias of the IF estimate using MBD 
have the same rates of change with respect to the window length h as those using 
WVD [9]. 

10.3.3.3 Examples of Quadratic TFDs Suitable for Multicomponent IF Esti- 
mation 

TFDs with time-only (or lag-independent) kernels constitute a subclass of the 
quadratic class of TFDs. These TFDs share the important properties of cross- 
terms suppression, high-resolution, and supporting amplitude estimation, making 
them well suited for multicomponent IF estimation. The modified B-distribution 
(a. k. a. hyperbolic T-distribution) was defined earlier in this section. Another 
example is the exponential T-distribution, which is defined in terms of its time-lag 
kernel as [8] 

G(t,r) = G v (t) = \fr\F ex p(~ ?7t 2 ) 

where ?y is a real parameter and \J ti/tt is a normalization factor. The resulting TFD 
used for multicomponent IF estimation is then given by eq. (10.3.4). 

10.3.4 An Adaptive Algorithm for Multicomponent IF Estimation 

Eq. (10.3.15) shows that the optimal window length using the MBD is a function of 
time and depends on the second derivative of the IF law /^(t); it decreases when 
the IF law f{t) has a high variation. Hence a time- varying window length is needed 
to optimize the estimation. The Stankovic-Katkovnik adaptive algorithm developed 
in [9] for monocomponent FM signals can be used since the IF estimation variance 
is a continuously decreasing function of h while its bias is continuously increasing, 
as shown in eqs. (10.3.13) and (10.3.14); see also Article 10.2. These conditions 
are necessary for bias-variance tradeoff such that the algorithm converges at the 
optimum window length that resolves this tradeoff. It is shown in [9] that, if h is 



(10.3.15) 
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small enough then the IF estimate fh(t) is inside the confidence interval D defined 
as follows 

D = [ fh(t ) - 2 Ky / var(A/ /l (t)), f h (t) + 2k\J va,r(A f h (t))\ (10.3.16) 

with Gaussian probability P(k), k being a parameter (usually 2); while for large h, 
fh(t) is outside D. Hence, if we consider an increasing sequence of window lengths 
{h r \r = 1 : TV} (TV being the number of samples) and calculate the MBD (and 
hence fh r (t)) for each h r , then all {D r } have at least one point in common (which 
is fh r (t)) ^ is sufficiently small. The first h r for which P r _i and D r have no 
point in common is considered optimal as it decides the bias-variance tradeoff. 

The estimates for the amplitude of the signal a and the variance of noise a 2 
[used in eq. (10.3.13) and implicitly in eq. (10.3.16)] were given in [9] as: 

i N i N 

a 2 + of = — | z(nAt) | 2 ; of = — ^ | z{nAt) - z((n - 1 ) At) | 2 (10.3.17) 

n—1 n—2 

where TV is the number of samples. For further details of this adaptive algorithm 
see [9]. 

For a multicomponent analytic signal of the form stated in eqs. (10.3.1) and 
(10.3.2), with {a m } constant, we can use the extension of the monoconrponent IF 
estimation algorithm in [9] for multicomponent signals as described in [2,3]. This 
algorithm tracks component maxima in the time-frequency plane and requires a 
threshold T p {t) so as to ignore the local maxima caused by the cross-terms and 
windowing. In fact, T p (t) is application and distribution dependent. 

The algorithm requires the knowledge of the confidence intervals D r m for each 
component, where r refers to the window length ( h r ) and m refers to the signal 
component. The calculation of -D r ,m depends on the estimation of the individual 
amplitudes a m of the components. Using the MBD, the actual amplitudes |a m | can 
be estimated as shown in [3]. Using |a m | 2 and a 2 to calculate var(A //(f)) [given by 
eq. (10.3.13) for MBD(f, /)], we can define the confidence intervals {D rtTn } for all 
components as in [2,3]. The IF f m (t) is contained in at least one of the confidence 
intervals { D r m } if h r is sufficiently small, and the optimal window length is the 
first h r (from the increasing sequence {h r \r = 1 : TV}) for which D r _ i >Tn and D r ^ m 
have no point in common. 

Example: We consider a three-component FM signal z(nAt) with amplitudes a± = 
0.5, 02 = 1, and 03 = 1.5 and non-linear IF laws: f± = 47+2.5 sinh - 1 (20(nAf— 0.2)), 
/2 = 30 + 2. 5sgn(40(nAt — 0.6)), and / 3 = 10 + 2 sin(10(?rAt — 0.7)), with SNR = 15 
dB, (3 = 0.1, k = 2, 0 < riAt < 1, and At = 1/128. Fig. 10.3.3 shows the result of 
the tracking adaptive algorithm for IF estimation of z(nAt) using the peaks of the 
MBD and the spectrogram. 

Fig. 10.3.4 shows the conventional peak IF estimation for the same signal using 
MBD and the spectrogram. Both TFDs fail to give a robust IF estimation at the 
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Fig. 10.3.3: Left: Adaptive IF estimation using the peak of the MBD for a three-component FM 
signal with total length N = 128, SNR = 15 dB, and At = 1/128. Dashed lines represent the true IF 
laws. Right: Adaptive IF estimation using the peak of the spectrogram for the same signal, assuming 
that component amplitudes are known. 
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Fig. 10.3.4: Left: IF estimation for the three-component FM signal as in Fig. 10.3.3 using the conven- 
tional (constant window) peak estimation. Left: MBD peak IF estimation. Right: Spectrogram peak 
IF estimation. (In addition to failure of IF estimation at instants of rapid frequency change, the spec- 
trogram has poor tracking performance. Note also that both methods cannot track the continuously 
varying frequency of the third component.) 



instants of rapid frequency change. In addition, the spectrogram has poor time- 
frequency resolution, both in adaptive and const ant- window IF estimation. 

10.3.5 Summary and Conclusions 

Concurrent IF estimation of the separate components of a multicomponent FM 
signal using TFD peaks location requires conditions on the selection of a suitable 
quadratic TFD. Required properties are: (1) high time-frequency resolution while 
suppressing cross-terms efficiently, (2) the TFD to enable direct amplitude estima- 
tion for the individual components, (3) the variance of the IF estimate using the 
TFD should be a continuously decreasing function of the lag window length while 
the bias is continuously increasing. Quadratic time-frequency distributions that 
satisfy these conditions were presented and discussed. A constant-window tracking 
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algorithm may not give a robust IF estimate if the IF changes rapidly with time 

due to the effect of the higher-order derivatives of the IF law. Hence an adaptive 

algorithm is used for robust multicomponent IF estimation. 
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